Growth of non-infinitesimal perturbations in turbulence 
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We discuss the effects of finite perturbations in fully developed turbulence by introducing a 
measure of the chaoticity degree associated to a given scale of the velocity field. This allows one 
to determine the predictability time for non-infinitesimal perturbations, generalizing the usual con- 
cept of maximum Lyapunov exponent. We also determine the scaling law for our indicator in the 
framework of the multifractal approach. We find that the scaling exponent is not sensitive to inter- 
mittency corrections, but is an invariant of the multifractal models. A numerical test of the results 
is performed in the shell model for the turbulent energy cascade. 
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The standard characterization of the chaotic behavior 
of a dynamical system is given by the maximum Lya- 
punov exponent A max , which measures the typical expo- 
nential rate of growth of an infinitesimal disturbance |Q . 
It is thus expected that the predictability time is pro- 
portional to A~ ax , the shortest characteristic time of the 
system. The underlying point is that the growth of a 
perturbation is well described by the linear equations for 
the tangent vector even if this cannot be literally true for 
non-infinitesimal perturbations. There exist indeed many 
situations where the Lyapunov analysis has no relevance 
for the predictability problem and it is necessary to in- 
troduce indicators which are able to capture the essential 
features of a chaotic system. For instance, when two or 
more characteristic time scales are present a direct iden- 
tification of the Lyapunov and predictability times leads 
to paradoxes as recently pointed out in Ref. . 

In this letter, we introduce a measure of the chaotic- 
ity degree related to the average doubling time that ex- 
tends the concept of Lyapunov exponent in the case of 
non-infinitesimal perturbations. Our indicator is a scale- 
dependent Lyapunov exponent which becomes particu- 
larly useful when there exists a hierarchy of characteristic 



times such as the eddy turn-over times in three dimen- 
sional fully developed turbulence || . 

In turbulent flows it is natural to argue that the max- 
imum Lyapunov exponent is roughly proportional the 
turnover time r of eddies of the size of the Kolmogorov 
length -q (the viscous cut-off) that is the shortest charac- 
teristic time pj. By dimensional analysis, the turnover 
time of an eddy of size £ is t(£) ~ where h is the 

scaling exponent of the velocity difference in the eddy 

vt=\v(x')-v(x)\~i h , l=\x'-x\. (1) 

The viscous cut-off vanishes as a power of the Reynolds 
number Re, i.e. rj ~ Re _1 ^ 1+/l - ) . These relations imply 
that the maximum Lyapunov exponent should scale as 

A max ~ 1 / 7-fa) ~ Re Q with a = (2) 

In the Kolmogorov K41 theory g, h = 1/3 for all 
space points x so that a = 1/2, as first pointed out by 
Ruelle Qj. However, the intermittency of energy dissi- 
pation leads to the existence of a spectrum of possible 
scaling exponents h affecting the value of a. In the mul- 
tifractal approach 0, the probability that the velocity 
difference on scale I scales as vi ~ £ is assumed to be 
P e (h) ~ £ 3 - D ( h \ This ansatz can be tested by measuring 
the scaling of the structure functions 
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In the K41 theory ||{| C P = p/3 while in the mul- 
tifractal scenario M £ p is a non-linear function of p 
given by the Legendre transform of the function D(h), 
C P = Tah\h[hp — D(h) + 3]. Moreover, as a conse- 
quence of multifractality there is a spectrum of viscous 
cut-offs M , since each h selects a different damping scale 
rj{h) ~ i?e~ 1 ^ 1+ ' 1 ' 1 , and hence a spectrum of turnover 
times Th(rj). To find the Lyapunov exponent, we have to 
integrate over the /i-distribution (|] 



p v {h)dh 



h - D ^ +2 dh~Re a . 



(4) 



In the limit Re — > oo, the integral can be estimated by the 
saddle point and gives a = maxj, [D(h) — 2 — h]/(l + h). 
By using the function D(h) obtained with the random 
beta model fit [p||lC|] one has a = 0.459... 

In the predictability problem, we are interested in 
defining the growth of an error on the velocity field. As 
usual we consider the Euclidean norm 



5v(t) 



d 3 x \v'(x,t) - v(x,t)f 



1/2 



(5) 



to introduce the notion of distance between two velocity 
fields v and v'. 

Then, the predictability time T p is the time necessary 
for an initial error Sv(0) = Sq to become larger than a 
given but arbitrary threshold value A: 



T„ = max [ 5v(t') < A for t' < t}. 
1 t>o 



(6) 



In a first approximation, neglecting the non-linear terms 
of the evolution equation for the error growth and assum- 
ing that both 5q and A are infinitesimal, one obtains 



T p ~ A-^ ln(A/|* | 



A: 



(7) 



In turbulence, such a relation would imply that the pre- 
dictability time decreases with the Reynolds number as 
Re _Q . This is contradictory with the quite intuitive ex- 
pectation that the predictability time should be roughly 
proportional to the turn-over time of the energy contain- 
ing eddies on the large scales, and so practically indepen- 
dent of the Reynolds number ptf^l . 

The paradox stems from assuming the validity of the 
Lyapunov analysis for perturbations Sv that are much 
larger than the typical velocity difference ~ rj/rij]) 
on the Kolmogorov length scale rj. In this case, the er- 
ror growth is non-exponential as it can be understood by 
simple heuristic arguments |TT|JT^ ]. The problem can be 
faced by generalizing the concept of maximum Lyapunov 
exponent to the case of non- infinitesimal perturbations. 



The generalization is particularly useful in systems with 
many characteristic time-scales. 

For this purpose, it is convenient to consider the time 
T r (5v) necessary for a perturbation to grow from Sv to 
rSv, for a generic r > 1. For r = 2 this is the doubling 
time of a perturbation, usually studied in atmospheric 
predictability experiments p^|Jl2] |. After performing an 
average over different realizations of the flow or, equiva- 
lently, a time average along a trajectory v(t) in the phase 
space, we introduce the scale-dependent Lyapunov expo- 
nent 



X(5v) 



1 



T r (Sv) 



In r. 
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Such a definition is consistent with the request of recov- 
ering the maximum Lyapunov exponent in the limit of 
infinitesimal error, since 



lim X(Sv) 

Sv^O 



(9) 



It is easy to estimate the scaling of X(Sv) when the per- 
turbation is in the inertial range iv, < Sv < vl, L 
being the size of the energy containing eddies. In this 
case, following the phenomenological ideas of Lorenz, 
the doubling predictability time of an error of magni- 
tude Sv can be identified with the turn-over time r(£) of 
an eddy with typical velocity difference vi ~ Sv. Since 
t(£) ~ £ h ~ 1 ~ v\ 1 ^ h , one has 



X(Sv) ~ Sv 



-13 



(3 = l/h - I. 



(10) 



Neglecting intermittency, i.e. 
value h = 1/3, gives (3 = 2. 



using the Kolmogorov 
In the dissipative range 



Sv < v rjl the error can be considered infinitesimal, imply- 
ing X{Sv) = A max . 

The intermittency of energy dissipation reflects the dy- 
namical intermittency of the chaoticity degree, so that 
our arguments based on dimensional analysis cannot be 
fully correct. In the framework of the multifractal ap- 
proach, our indicator scales as 



X(Sv) 



dh Sv [3-n{h)Vh 8v i-i/h 



(11) 



where we have used arguments similar to those leading to 
([|) and the scaling factor £ ~ Sv 1 / 11 . ^From the inequal- 
ity D{h) < 3h + 2, which is the analogous for turbulence 
of the standard inequality f(a) < a in multifractals, we 
have 



2 + h-D(h) 



> -2 for all h. 



(12) 



Equality holds for h = h 3 , the exponent that real- 
izes the minimum in the Legendre transform for the 
exponent of the third-order structure function £3 = 
min^ [3/i + 3 — D(h)] = 1. Therefore a saddle point esti- 
mation of (pd|) gives 
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j3 — min 

h 



2 + h- D(h) 



= -2. 



(13) 



An important consequence of multifractality follows 
from the existence of a spectrum of dissipative cut-offs 
r](h) which reduces the effective inertial range where the 
scaling of X(Sv) holds [[s|Jl3|] . To be more specific, the mul- 
tifractal approach leads to rj(h) ~ Re~ 1 ^ 1+,l - ) , and the 
integral (|ll| ) has to be performed for h(5v) < h < h max , 
where h(5v) is given by 



h(8v) 
6v~Re 1 + h (Sv) 



(14) 



As a consequence, the scaling X(5v) ~ 8v P holds only for 
5v > Re~ ,l3 ^ 1+h3 \ i.e., the inertial range is reduced by 
intermittency. In the range Re -1 / 4 < 8v < R e ~' l3 /( 1 +' 13 ) 
we expect a non-trivial shape of X(Sv) depending on 
D(h). 

In order to test our results we have numerically stud- 
ied the GOY shell model (l4],|l5| for the energy cas- 
cade in fully developed turbulence. The model is an 
approximation of the Navicr-Stokes equations obtained 
by dividing the Fourier space into shells of wavenumbers 
k n < \k\ < k n+ i. A complex scalar u n is associated with 
the n th shell individuated by k n = k 2 n . It represents 
the velocity difference over a length scale £ ~ fc" 1 . Since 
the energy cascade in turbulence is believed to be local 
in the A;-space with an exponentially decreasing interac- 
tion among shells, it is reasonable to consider only the 
interactions of a shell with its nearest and next-nearest 
neighbors. The Navier-Stokes equations are then approx- 
imated by a set of ordinary differential equations: 



d_ 

dt 



u n = g n - vk n u n + f5 nA 



(15) 



9n — i a nk n u n j r iU n j r 2 + ib n k n —\u n _}U n j r Y + 

ic„fc„_ 2 <_2<-l ( 16 ) 

with b\ = &at = ci = C2 = djv-i = ayv = 0. The coeffi- 
cients of the nonlinear term obey a n + + c„+2 = 
to ensure energy conservation for / = v = 0. With the 
standard choice for three dimensional turbulence a n = 1 , 
b n = —1/2 and c„ = —1/2, there is a second conserved 
quantity (— l) n k n \u n \ 2 , which in the shell model plays 
the role of helicity [[l6| . 

The shell model exhibits non-linear exponents £ p for 
the structure functions as found in experimental 

data M ■ We have determined the scale dependent Lya- 
punov exponent by a numerical integration of the GOY 
model starting form two different initial conditions. The 
distance between the two velocity fields is Euclidean, 
i.e. Su = (%2 n \u n — u^J 2 ) 1 / 2 and we have computed a 
time average over the trajectory {u n (t)} of the quantity 
TriSu)- 1 with r = 2 1 / 2 . 

Figure [l] shows the scaling of (\/T r {8v)) as a function 
of 5v in the GOY model with N = 27 shells and viscosity 



v = 10 9 . For comparison we have also computed the 
eddy turn-over times 



k n {\ u i 



,2x1/2 



(17) 



We thus provide a clear evidence that there is a large 
range of small scales where X(Su) — A max and r n ~ , 
as a consequence of the reduction of the inertial range 
for the scale dependent Lyapunov exponent. Note that 
in the GOY model A max ~ 10 _2 r(7y)~ 1 , though the de- 
pendence of the two quantities on the Reynolds number 
is the same. 

In Figure ||, we have plotted the rescaled quantity 
X(Su)/Re a versus Su/Re~ 7 using the Kolmogorov val- 
ues h = 1/3, a = (1 - h)/(l + h) = 1/2 and 7 = 
h/(l + h) = 1/4. The collapse of the data obtained at dif- 
ferent Reynolds numbers is a pretty confirmation of our 
scaling arguments. The reduction of the inertial range 
for the scale dependent Lyapunov exponent reveals itself 
in the small range where X(5v) ~ Sv~ 2 holds. 

We conclude by noting that our scale-dependent Lya- 
punov exponent X(Sv) has some similarity with the con- 
cept of the e-entropy recently discussed by Gaspard and 
Wang [p7|-^9| for the treatment of experimental data. 
We stress that in chaotic systems the maximum Lya- 
punov exponent is often more relevant and more easily 
computed than the Kolmogorov- Sinai entropy. There- 
fore we believe that our X(5v) will often be more relevant 
and more easily computable than the e-entropy. More- 
over, since we use the evolution law and not experimental 
data, in our case there are no particular limitations as to 
the number of degrees of freedom involved. 

In conclusion, when the perturbations are non- 
infinitesimal it is necessary to extend the definition of 
Lyapunov exponent to make it physically consistent. The 
generalization proposed in this letter is particular useful 
when many characteristic time scales are present. Our re- 
sult allows one to get a quantitative control of the growth 
of perturbations which are non-infinitcsimal, looking at 
the average of the inverse doubling time. By this def- 
inition one has the two advantages of maintaining the 
link with the forecast limitation of a system and of re- 
covering the maximum Lyapunov exponent in the limit 
of infinitesimal perturbations. The scale-dependent Lya- 
punov exponent thus is an important tool of investigation 
of highly dimensional dynamical systems and, far from 
being limited to the predictability problem of turbulent 
flows in geophysics [[ll 12 20]] , it can assume a great rel- 
evance in the characterization of very different chaotic 
phenomena. 
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FIG. 1. (1/T r (8v)} (diamond) as a function of 8v for the 
COY model with N = 27, k = 0.05, / = (1 + i) x 0.005 and 
v — 10 . The crosses are the inverse of the eddy turn-over 
times t _1 (5w) = k n (\u n \ 2 ) 1/2 versus Sv = (|u„| 2 ) 1/2 . The 
straight line has slope —2. 



FIG. 2. In [(1/T r (,5«))/Re 1/2 ] versus In [Sv/Re' 1 ^] at dif- 
ferent Reynolds numbers Re = v~ x . The results are ob- 
tained in the COY model for k = 0.05, / = (1 + i) x 0.005 
and: (diamond) N = 24 and v = 10~ 8 ; (plus) N = 27 and 
v = 10~ 9 ; (square) N = 32 and v = lO" 10 ; (cross) N = 35 
and v — 10~ n . The straight line has slope —2. 
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